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Abstract. We use cosmic microwave background (CMB) data from the 9-year 
WMAP release to derive constraints on monodromy inflation, which is characterized 
by a linear infiaton potential with a periodic modulation. We identify two possible 
periodic modulations that improve the fit, by A^^ ~ 10 and 20 respectively. However, 
standard Bayesian model selection criteria assign roughly equal odds to the modulated 
potential and the unmodulated case. A modulated inflationary potential can gen- 
erate substantial primordial non-Gaussianity with a specific and characteristic form. 
For the best-fit parameters to the WMAP angular power spectrum, the correspond- 
ing non-Gaussianity might be detectable in upcoming CMB data, allowing nontrivial 
consistency checks on the predictions of a modulated inflationary potential. 
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1 Introduction 

Monodromy inflation is a theoretically well-motivated scenario for implementing a 
period of nearly exponential expansion in the very early universe [1-3] . Derived within 
string theory, monodromy models "unwrap" angular directions in fleld space, and have 
monomial potentials over super-Planckian distances. The underlying periodicity of the 
fleld may be imprinted on the potential in the form of a small oscillatory modulation. 
The evolving inflaton is often described as a clock that coordinates the inflationary 
expansion over super-Hubble scales. The periodicity of the compact direction adds 
another timescale to the problem - so that the metaphorical inflaton clock now has 
both an hour hand and a second hand. 

Phenomenologically, monodromy inflation has a rich set of observable signatures. 
In particular, the Lyth bound [4] implies that the super-Planckian fleld ranges are 
associated with a detectable primordial gravitational wave background. The modu- 
lated potential naturally leads to an oscillatory power spectrum, which is a non-trivial 
signature of a monodromy-driven inflationary era. Moreover, monodromy inflation pro- 
vides a stringy implementation of resonant non-Gaussianity [5, 6], generating a 3-point 
function whose properties would be strongly correlated with the power spectrum. If 
couplings to other degrees of freedom are small, monodromy potentials have the ability 
to generate oscillons - long lived, localized, pseudo-stable excitations of a scalar fleld - 
in the post- inflationary universe [7]. Lastly, couplings of the inflaton to a gauge fleld, 
when present, may also lead to interesting phenomenology [8]. 

We determine the observational constraints from the cosmic microwave back- 
ground [CMB] on the canonical model of monodromy inflation, using the 9-year Wilkin- 
son Microwave Background Anisotropy Probe [WMAP] dataset [9, 10]. We estimate 
the likely amplitude of the 3-point signal associated with the parameter values that best 
flt the angular power spectrum data. The speciflc form of this signal is not well cap- 
tured by standard measures of non-Gaussianity, but can be recovered from CMB data 
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with a suitable estimator, such as a modification of the modal expansion method [11] 
(see, e.g., Ref. [12]). We show that information from the 3-point function could tighten 
constraints on any possible modulation. 

Our analysis uses the ModeCode package^ [13-15], an extension to CoSMoMC 
[16] developed by two of us [RE and HP] and other collaborators. ModeCode solves 
the inflationary perturbation equations numerically using an algorithm originally devel- 
oped in Ref. [17], avoiding any analytical approximations beyond those inherent in the 
underlying cosmological perturbation theory. This approach allows us to directly esti- 
mate the free parameters in the inflationary potential, independently of any empirical 
characterizations of the primordial power spectrum. Separately, ModeCode replaces 
the usual Markov Chain sampler in CoSMoMC with the nested sampler MultiNest 
[18, 19]. This copes naturally with multimodal likelihood surfaces, which can arise with 
modulated potentials or "feature" models, and simultaneously computes the Bayesian 
evidence (e.g., Ref. [20]). 

Calculating the Bayesian evidence addresses the model selection problem, deter- 
mining whether the improved fit to the data that follows from inclusion of the modu- 
lations - or any other parameters we might add - justifies the additional complexity of 
the resulting model. We compare the modulated potential to the standard monodromy 
potential without modulations, finding that the usual model selection criteria assign 
slightly better odds to the modulated model, despite the fact that the former has three 
extra parameters. 

This paper is organized as follows. In Section 2 we summarize the properties of 
monodromy inflation and specify the priors for the parameters in the potential. In 
Section 3 we present the parameter estimates and Bayesian evidence values. Section 4 
discusses the expected 3-point statistics that are permitted by the constraints derived 
for the potential, and we discuss our results in Section 5. 

2 Monodromy Inflation, Perturbations and Methodology 

2.1 Potential and Model 

We work with the monodromy potential 



Here and / are energy scales, while 6 is a dimensionless parameter and we include 
a constant phase as in Ref. [3]. For convenience, we have introduced an additional 
dimensionless parameter c and choose it so that the potential vanishes at = 0. This 
potential describes the inflationary phase but is unbounded below, and thus does not 
describe the reheating physics. We assume that (p is initially large and positive. 



^ ModeCode can be downloaded from http://zuscrver2.star.ucl.ac.uk/^hiranya/ModeCode/ and 
applies a patch to CosmoMC. For further details and instructions consult the documentation supplied 
with the code. 




(2.1) 
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This potential is simplified relative to the "full" monodromy potential in Ref. [2], 
which has a quadratic minimum at the origin. The version here does not significantly 
modify the potential in the region relevant to the cosmological perturbations, but does 
change the field value at which inflation ends. 

The equations of motion obeyed by the field, scale factor and perturbations as 
implemented in ModeCode are described in Refs. [13-15]. We will need to refer to 
the number of e-folds of inflation 



N = / ^ ^ iV ~ ;rv7^ . (2.2) 

The approximate equality is a reliable estimate for moderate values of b. For = 55, (p 
is a little larger than lOMpi and between (say) = 55 and = 45, 6(j) ~ Mpi. Setting 
aside the physical constraints on / described below, if / ^ 10~^Mpi the portion of 
the potential which generates perturbations on astrophysical scales will resemble an 
unmodulated potential with a "running" spectral index (see also Ref. [21]) while if 
/ ^ Mpi the modulations will cause the power spectrum to vary on a scale that is 
very short compared to all the other variables in the problem. 

2.2 ModeCode, Power Spectrum and Evidence 

The monodromy power spectrum can be computed via an approximate but accurate 
analytic calculation [3, 6] or numerically via the routines in ModeCode, which easily 
accommodates new potentials. We worked with the numerical calculation as it remains 
valid even if the model parameters take extreme values. In order to accurately compute 
the power spectrum in the small / limit, the initial conditions for the mode evolution 
are set further inside the horizon than is required for other models, and we significantly 
boosted the sampling of the primordial power spectrum P{k) and the accuracy settings 
within CAMB, relative to their default values. 

ModeCode incorporates the MultiNest sampler and computes the evidence 
E, the integrated likelihood £ over the parameter volume {ai, ■ ■ ■ , a^/ }, 

E = j da^' P{ai)C{ai) , (2.3) 

weighted by the prior P(aj), normalized so that j da^^ P{ai) = 1. For uniform priors 
P{ai) is a constant (or zero), and the weighting is the inverse parameter volume, or 

E=-^ [ da^'Cia,) . (2.4) 
VoIm j 

Evidence allows comparisons between two models, rather than an absolute measure of 
quality. Our comparison model is inflation with a purely linear potential V{(j)) = fi^(j). 
The computed evidence depends on the parameter ranges allowed by the priors, as 
discussed in detail in Ref. [14]. 
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Inflation 


Mass scale 


-3.615 < logio(/i/Mpi) < -3.015 


Axion decay constant 


-3.4 < logio(//Mpi) < -2.0 


Oscillation amplitude 


< 6 < 0.9 


Phase 


— 7T < Ip < 71 


Matching 


e-foldings 


N = 55 


Astrophysics 


Bar yon fraction 


0.0218859 < Qhh'^ < 0.02378859 


Dark matter 


^dmh^ = 0.1145 


Reionization 


r = 0.0874 


Projected acoustic scale 


e = 1.040 


Sunyaev-Zel'dovich Amplitude 


Asz = 0.10078 



Table 1. Parameter priors used in the analysis. For linear inflation, 6 = and the /, ■0 
parameters are unused. We assume that ^Ik is zero in all cases, and parameters are drawn 
from a uniform prior. The pivot scale is k^, = 0.05 Mpc^^. The non-primordial parameters 
except for ilb^^ ai'e fixed at their maximum likelihood values for the unmodulated case. 

2.3 Model Specification and Priors 

The energy scale of thermalization, and thus the "average" equation of state in the 
post-inflationary universe, is largely unconstrained - there is no hard evidence that 
the universe is thermalized above MeV scales at which neutrino freeze-out occurs. The 
mapping between present day scales and the value of at which the corresponding 
mode leaves the horizon and the value of at horizon crossing is thus poorly constrained 
[14, 22]. Consequently, a (very) small change in the post-inflationary equation of state 
is degenerate with the phase variable in the potential; therefore, in what follows, we 
will flx the number of e-folds at which the pivot scale leaves the horizon to 55. 

A preliminary investigation of this system varying both the standard cosmological 
parameters and the free parameters in the potential showed that the posterior on 
logio (//^Pi) has a peak near / ~ Mpi, approaches zero near log^^g (//^pi) ~ —1.5 
and rises substantially at smaller values of /. Given that values of / > 10~^Mpi are 
inconsistent with the assumptions under which the monodromy potential was derived, 
it is clear that we only need to concern ourselves with a "fast" modulation of the 
power spectrum in the small / limit. Consequently, we can perform our analysis after 
fixing the usual cosmological parameters, except for Q^^h"^ which potentially exhibits a 
slight degeneracy with / (see below). The specific values we use to fix the remaining 
parameters came from the an initial run performed without modulations, and the priors 
are presented in Table 1. 

Physically, only weak constraints can be placed on the free parameters in equa- 
tion (2.1). The value of /i sets the amplitude of the primordial perturbation spectrum, 
but is not predicted by fundamental theory. As described in Ref. [14], ModeCode ex- 
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Figure 1. Marginalized posterior distributions for inflationary parameters. 



eludes parameters for which the power spectrum amplitude is incompatible with basic 
structure formation requirements, and this stipulation dominates the permitted range 
of n when it is specified with a logarithmic prior. The axion decay constant / is also an 
unknown mass scale with a large possible range. On very general grounds we expect 
/ < Mpi [23, 24], and the approximations underlying the derivation of the modulated 
potential become unreliable if / ^ 10~^Mpi [2, 3]. For a small range / > 10~^Mpi 
there is a degeneracy between / and f2b [3]. To avoid such potential degeneracies 
with the standard cosmological parameters, we work with / < 10~^Mpi, and we vary 
f2b^^ in order to eliminate any effects arising from this degeneracy. We can then self- 
consistently fix the remaining parameters to their best-fit values. The validity of the 
approximations made in Refs. [2, 3] also imposes a lower bound on /, and we use 
/ > 4 X 10~^Mpi. (A similar bound can also be derived by requiring weak coupling of 
the effective theory [25].) In order to avoid potential numerical singularities, we work 
with the limit b < 0.9. 



3 Parameter Estimates and Evidence 

We base our constraints on the WMAP 9- year dataset [9, 10]. Constraints on the 
monodromy model from WMAP7 were derived in Refs. [26, 27], and other studies of 
oscillatory features in the CMB include Refs. [28-31]. The Atacama Cosmology Tele- 
scope [ACT] [32] and South Pole Telescope [SPT] [33] datasets contain information 
on the power spectrum at angular scales unresolved by WMAP. However, the ACT 
and SPT likelihoods are binned into bandpowers, limiting their sensitivity to short 
wavelength modulations. Moreover, the two datasets lead to apparently different con- 
clusions about the form of the power spectrum in the smooth limit [34] . Consequently, 
we do not include either of them in our analysis. Modulations in the primordial power 
spectrum also lead to modulations in the galaxy power spectrum measured from large 
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Figure 2. Reference angular power spectrum with 6 = and difference for the best fit 
spectrum with logj^g (//-^Pi) ^ —3.38 in blue and with log;^o (f/Mpy) —3.18 in orange. 

scale structure surveys. Currently these bounds are weaker^ and we focus on the 
WMAP data. 

Figure 1 shows posterior distributions for the free parameters / and b. Two spikes 
in the posterior for log^Q (//Mpi) are visible, and the improvement in the effective or 
— 21og£ associated with the lefthand peak (log^Q (//Mpi) ^ —3.38) is approximately 
Axlff = 19.4 over the unmodulated potential, and Axis — 18.9 with respect to ACDM, 
which has best-fit xls ~ 7558.0 for WMAP9. The posterior for b contains a peak away 
from 6 = associated with the two likelihood peaks; however, the data do not exclude 
6 = 0, and the marginalized posterior probability exhibits another rise there. We do 
not plot the posterior for the phase ip, as it is almost uniform, essentially reproducing 
the prior. This is consistent with the peaks seen in log^g (//Mpi), as the marginalized 



^See however Ref. [35] for a Euclid forecast. 
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posterior for the phase is not dominated by the narrow peaks in log^Q (//Mpi). 

The righthand peak at log^Q (//Mpi) ^ —3.18 leads to an improvement Axis — 
11.7 with respect to the unmodulated potential. It is consistent with the best-fit point 
in Ref. [3], where it had an improvement in the Ax^g of 11. Ref. [3] also finds a spike 
in the WMAP5 likelihood which is consistent with the higher peak seen here, in both 
frequency and amplitude. Ref. [3] also sees several smaller peaks, and these are less 
prominent in the new analysis. The two most prominent examples thus "survive", 
albeit after trading roles. The analysis in Ref. [3] made use of the analytic expression 
for the primordial power spectrum, which provides a non-trivial consistency check for 
ModeCode. When b is increased at fixed fi and / ^ Mp\, the averaged value of 
P{k) increases. Consequently, at large values of b, the best fit spectrum has a slightly 
smaller value of fi than the best fit in the 6 = limit. 

A number of other "feature" models are known to yield an improvement in Ax^ of 
order 10 (e.g. a step in the inflaton potential, yielding high-£ oscillations in the power 
spectrum [36]) and even if one of them were to be genuine, they cannot all be genuine. 
However, the enhancement in the fit at the logj^g (//^pi) ~ —3.38 peak is relatively 
substantial. We have performed preliminary simulations to determine the significance 
of such an improvement in a fit to data without a signal, which indicate that the 
improvement is less significant than it may seem at first sight, in accordance with the 
results from the model selection analysis. Further, the approximations in the WMAP 
likelihood code have not necessarily been tested in the regime of very high frequency 
oscillations in which our best fit parameters lie. The ACg for the best-fit parameters 
associated with the higher peak in the posterior distribution of logiQ(//Mpi) are shown 
in Figure 2. In Figure 3 we show the contributions to Axis ^ function of multipole, 
for the higher likelihood peak. We see that the reduction in comes from the full 
multipole range, rather than being concentrated in a narrow region, which is consistent 
with a modulated model. 

The presence of sharp features in the marginalized posteriors taxes the usual 
software tools used in cosmological MCMC analyses. Fortunately, we can make use 
of the MultiNest sampler [18], which copes with multimodal and rapidly varying 
likelihood functions. Plotting binned 2D posteriors leads to messy results, and we 
instead show the "heat map" of the posterior in Figure 4, mapped to a qualitative 
intensity scale. 

We compute the evidence ratio between the pure, unmodulated monodromy po- 
tential and the modulated case using MultiNest, finding Alog-E ~ —0.6 in favor of 
the latter. Thus, the betting odds favor the modulated monodromy model by about 
2:1 despite its three extra parameters. However, these odds are far from decisive^, and 
the WMAP9 data cannot differentiate between these competing models. 
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Figure 3. The difference in Ax^g given by the WMAP9 MASTER TTTT likehhood for 



the reference and best-fit modulated models as a function of 
labels. 



binned as indicated by the 



4 The Three Point Function 



Up until this point much of our analysis is functionally equivalent to working with a 
modulated primordial power spectrum P{k) (e.g. Ref. [28, 37, 38]). However, mon- 
odromy inflation is a stringy realization of the resonant non-Gaussianity scenario, first 
described in Ref. [5], and makes unambiguous predictions for both the 2-point func- 
tion and the 3-point function. Consequently, given our posterior we can also compute 
the expected amplitude of the 3-point function. The "shape" of this non-Gaussianity 
does not map directly onto any of the forms with known constraints, so we cannot use 
the CMB 3-point function in the parameter estimation process. However, we can ask 
whether this information would tighten constraints derived from the 2-point function 
or, more optimistically, confirm that a peak seen in the posterior of log^^g (//-^pi) is 
consistent with an initial period of monodromy inflation. 

For this model the approximate amplitude of the 3-point function is [6] 



3V27r6 



NL 



(4.1) 



^Betting odds of 150:1 are considered decisive in this context. 
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WMAP9 Monodromy 




-3.35 -3.30 -3.25 -3.20 -3.15 -3.10 -3.05 

logio(f/Mpi) 

Figure 4. Marginalized 2D-joint posterior for b and log^^Q (//Mpi) - the intensity is a 
function of the posterior, and shows that the peaks are narrow and well-defined. Note that 
the color scale is qualitative. 

where 0^ is the field value at which the perturbations are generated, which is roughly 
given our priors. We plot the inferred posterior for /nl in Figure 5 - this 
is peaked strongly near zero, but integrating the tail we find that -P(/nl > 100) ~ 0.75. 

To our knowledge, there are no published constraints on /nl for resonant non- 
Gaussianity, or even estimates of the likely constraining power of datasets such as 
Planck^ for this range of frequencies. However, a simple estimate presented in Ref. [25] 
indicates that the signal-to-noise ratio for a measurement of the 3-point function is 
smaller but comparable to the one in the 2-point function in the regime of interest. 
We plot the predicted value of /nl as a function of b and log^g (//^pi) i^ Figure 6, 
and overlay the best fit points for the two peaks seen in the posterior for log^^g (//^pi)- 
Resonant non-Gaussianity will be more challenging to constrain than some of the 
simpler forms. However, future data may constrain the power spectrum and /nl with 
sufficient accuracy to test whether the peaks seen in the posterior for log^^Q (//Mpi) are 



^http:/ /www.esa.int/Planck 
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Figure 5. Posterior distribution for /nl- 



associated with a modulated inflaton potential. Further motivation for a measurement 
of the 3-point function comes from models presented in Ref. [39] that can lead to 
resonant non-Gaussianity without a corresponding signal in the power spectrum. 



5 Discussion 

We have used ModeCode [13-15] to derive constraints on monodromy inflation based 
on the WMAP 9- year dataset. We find two prominent peaks in the marginalized 
distribution for log;^o (//^pi); the axion decay constant. 

The improvement in the at the peaks is ~ — 12 and ~ —20 with respect to the 
unmodulated potential, and slightly worse with respect to ACDM. The concordance 
cosmological model is a very good fit to the WMAP9 data, with per degree of 
freedom 1200/1168 for £ = 33 - 1200, i.e., ~ 1.03 [9]. Thus, a model which better 
describes the data has access to an improvement of Ax^g ~ 30 on average, before it 
begins over-fitting the data. The higher likelihood peak in the monodromy case gives 
per degree of freedom 1183/1168 in the same multipole range, i.e., ~ 1.01. 

The monodromy model contains a single modulation, so within the context of our 
prior, at least one of these peaks is a spurious "fit to noise" , rather than evidence for a 
genuine modulation of the primordial power spectrum. The improvement in fit at the 
higher peak in the posterior for log^g / is noticeable, but not compelling. Moreover, 
even if the modulation was present in the processed data, additional work would be 
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Figure 6. Contours show the predicted amphtude of the resonant non-Gaussianity signal 
associated with the the modulation parameters / and b. The two points show the best fit 
values associated with the peaks seen in the marginalized distribution for log^o (//Afpi)- 

required to demonstrate that this was not an artifact of the data analysis.^ 

In future CMB datasets, the polarization is likely to be tightly constrained. The 
scalar power spectrum dominates the (EE) and (TE) correlations and so any features 
in the underlying spectrum will produce a similarly modulated signal in the ii^-mode 
polarization, allowing much tighter constraints on any putative signal. Moreover, the 
polarization transfer function is narrower than the temperature transfer function, so 
features in the primordial spectrum are more clearly resolved in (EE) than in (TT) [41]. 
Beyond improved measurements of the power spectrum, we have also considered the 
amplitude of the 3-point function associated with monodromy potential, demonstrating 



^The first Planck cosmology data release took place after this work was first posted to the ArXiv. 
An initial search was made for a modulated power spectrum [40], but that scenario does not match 
the monodromy model considered here. 



that this is substantial for parameter values where the likelihood is peaked. If either 
of the candidate modulations seen here was confirmed within a higher-quality dataset, 
the predictions for the 3-point function are likely to be testable. This suggests that 
after detecting a modulated power spectrum - which would be an important result in 
its own right - it may be possible to determine whether this modulation was associated 
with a monodromy-style inflaton potential. Finally, monodromy inflation predicts the 
existence of a detectable background of gravitational waves, providing a further test 
of this scenario. 
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